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Abstract 

We perform a first principle calculation on the electronic properties of carbon doped graphitic 
boron nitride (graphitic-BN) . It was found that carbon substitution for either boron or nitrogen 
atom in graphitic-UN can induce spontaneous magnetization. Calculations based on density func- 
tional theory with the local spin density approximation on the electronic band structure revealed 
a spin polarized, dispersionless band near the Fermi energy. Spin density contours showed that 
the magnetization density originates from the carbon atom. The magnetization can be attributed 
to the carbon 2p electron. Charge density distribution shows that the carbon atom forms covalent 
bonds with its three nearest neighbourhood. The spontaneous magnetization survives the curva- 
ture effect in BN nanotubes, suggesting the possibility of molecular magnets made from graphitic- 
BN. Compared to other theoretical models of light-element or metal-free magnetic materials, the 
carbon-doped BN are more experimentally accessible and can be potentially useful. 
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Tranditional magnetic metal based materials have a Curie temperature not over 1,400 K. 
Thus there has been a long standing search for novel high temperature magnets Q|. Potential 
magnetism from metal-free system was probably firstly proposed by Fujita et al. Q]. they 
suggested the possible magnetic structure in the nano-scaled zigzag graphite ribbon resulting 
from the effect of localized edge states. The weak ferromagnetism in all-carbon system were 
observed by Makarova and Esquinazi et al. 4j respectively. These studies stimulated 
the interest to explore the possible mechanism of the unexpected magnetism in matal-free 
system. Lehtinen et al. fj] performed an ab initio local-spin-density approximation (LSDA) 
calculation to study the properties of a carbon adatom on a graphite sheet and found that 
this defect is magnetic with a moment up to 0.5 hb- This was important in understanding the 
magnetism experimentally observed by Esquinazi J] in the graphite system. More recently, 
Ma et al. [f| studied the magnetic properties of vacuncies in graphite and carbon nanotubes. 
They found that the vacancy was spin-polarized with a magnetic moment of 1.0 hb m the 
graphite sheets and also introduce magnetism in the carbon nanotubes, depending on the 
chirality of the nanotubes and the structural configuration with respective to the tube axis. 

The unexpected magnetism in metal-free system also ignited ambitions to search for 
metal-free magnetic materials in view of their potential applications as high temperature 
magnetis since metal magnets lose their ferromagnetism at high temperatures. Based on 
the results of Fjuita em et al. ^, Kusakabe and Maruyama \l\ further predicted that the 
nanoscaled zigzag-edged graphite ribbon could also have finite magnetization when each 
carbon at one edge is hydrogenated by two hydrogen atoms while each carbon atom at the 
other edge is hydrogenated with a single hydrogen atom. Another recently proposed model 
by Choi and coworkers j^] is the heterostructured C-BN nanotubes. They calculated the 
electronic structure of the (9,0) Ci(BN)i and C2(BN)2 nanotubes using density functional 
theory and found the occurrence of magnetism at the zigzag boundary connecting carbon 
and boron nitride segments of tubes. However, in view of their structural configurations, 
experimental synthesis was a great challenge. For Kusakabe's model, fully hydrogenating 
one edge while halfly hydrogenating the other side uniformly could not be done in the same 
experimental condition, what's more, in the half hydrogenation, it was hard to prevent 
two hydrogen atoms hydrogenating a single carbon atom which would destroy the magnetic 
ordering. As for the model of Choi js| , it required a BN nanotube segment to replace carbon 
compartment along the axis direction to form a C-BN nanotube superlattice which is very 
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hard in an enviroment where the three elements coexist since they could form covalent 
compound with each other. 

However, all these theoretical calculation had throught some hints about the nature 
of the magnetsim in the carbon-based system. That is, the magnetism may result from 
undercoordinated carbon atoms or unpaired carbon 2p electrons |7|, |8( . In this report 
we explore the electronic properties of carbon-doped graphitic-BN . The reason lies in the 
follwing reasons. First, graphtic-BN has a melting point over 1,700 K, which make it 
suitable for high temperature applications; Second , qr aphitic-BN can form different nano- 
scaled structures such as nanotubes Q. fune^Q and M1 en teS Q. The BN- b a 8 ed 
magnetic nano-structures can serve as molecular magnets; Third, most important of all, in 
the hexgonal network, a carbon atom may form three bonds with either boron or nitrigen 
atoms, leaving one electron unpaired (the interactions between layers are only weak Van der 
Waals, so the effect can be ignored) , making it a possibility of spontaneous magnetization 
in such a carbon containing boron nitride system. 

Our calculation is base on the spin-polarized density functional theory (DFT). The gener- 
alized gradient approximation (GGA) of Perdew and Wang is adaptted for the exchange- 
correlation potential. The projector augmented wave (PAW) potentials have been used to 
represent the electron-ion interactions. All 2s and 2p electrons are considered as valence 
electrons. A kinetic energy cutoff of 520 eV is used to ensure a convergence better than 1 
meV for total energies. A (3x3 xl) graphitic-BN sheet seperated by a 12 A vacuum under 
periodical bourdary conditions (PBC)is used for calculations. A Gamma-centered (4x4x1) 
K-mesh was used to sample the irreversible Brillouin Zone (IBZ). The conjugate-gradient al- 
gorithm jl^] is applied to relax the ions into their intantaneous ground states. The structural 
parameters have been adjusted so that the external pressure is less than 1 kBar along each 
lattice direction. We first calculate the structural and electronic properites of graphitic-BN 
and compare them to previous theoretical calculations and experimental studies before the 
calculations on the carbon doped supercell. All the calculation was done using the plane 
wave basis VASP package [l^ . 

We first calculated the lattice parameters and band structure of graphtic-BN sheet. Our 
calculation reached a lattice constant of 2.52 A and an indirect band gap of 4.61 eV. This 
reproduced previous calculation by Xu Q]. The validified the application of the projector 
augmented wave (PAW) potentials. The formation energy Ef orm is defined as 
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Eform = -E'toi(doped)-i?tof(perfect)-(£' y i(isolated)-i?B(isolated) 

where £U (isolated) and .Eg (isolated) is the energy of doping atom and the substituted 
atom at isolated states respectively. The calculated Ef orm is 1.95 eV for carbon substitution 
upon boron and 3.46 eV for carbon substition upon nitrigen. Compared to the energy 
required to dope a silicon atom substitution upon a carbon atom in Cqq which is 5.77 
eV and the experimental observation of CsgSii ^|, the energy required for a carbon atom 
to substitute boron or nitrogen atom in graphtic-BN is still within experimental access. 

Fig. 1 (a) is the band structure for the graphitic-BN with a boron atom substituted by 
carbon atom. We can see that the carbon substitution caused a lift up of the Fermi level of 
the system and the occurance of a flat band near Fermi level. This is because the carbon 
atom acts as a n-type dopant since a carbon atom consists 4 valence electrons. All bands 
with energies lower thatn 4 eV are fully occupied and thus have no contribution to the 
spin-polarization. However, the flat band near Fermi level is split into two branches. The 
spin-up branch is occupied and the spin-down branch is left empty, leading to a spontaneous 
polarization in the doped system with a net magnetic moment of around 1 /i^. The same 
effect also survive when a nitrogen atom is substituted by carbon atom (Fig. 1 (b)). The 
only difference is that the Fermi level is push down ( in this case, the carbon atom acts 
as p-type dopant because it have less one electron than nitrogen atom) and the occupied 
branch merges into those fully occupied bands. 

To find out the origin of the magnetism, we plot the magnetization density (spin-up minus 
spin-down density) contours through the sheet for both system in Fig. 2. As we can see that 
the magnetization density originates from where the carbon atom resides. No magnetization 
density is found at other atomic sites. This suggests that the spin-polarization results from 
the unpaired electron in the carbon atom. In Fig. 3, we plot the spin density distribution 
contour perpendicular to the sheet plane. In both system the spin density showed a p-orbital 
character. These results are reasonable as in the hexagonal network, a single carbon atom 
has three nearest neighbourhood to form bonds with, leaving one electron unpaired. It is 
obvious that in both case, the magnetization is from one of the carbon 2p electrons. 

One question in concern is that whether the carbon dopant can form bonds with its 
three nearest neighbourhood and maintain the fine structure. There has been successful 
omental synthesis of the graphitic BC3 _h epitaxial sheetfi We plot the 
charge density distribution contour through the carbon doped BN sheet in both cases shown 
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as Fig. 4 (a) and (b). In both systems we can see that carbon dopant forms bonds with 
nearest neighbourhood. When carbon is surrounded by nitrogen atoms (a), the electron 
clouds are localized around the connecting lines between carbon and nitrogen atoms. Due 
to the difference of electronegativity, carbon looses part of it valence electrons. In the case 
of substitution for nitrogen (b), carbon atom attracts electrons from boron atoms. In both 
case, the nature of the bonds is covalent and no structural distortion has been observed 
opun geometry optimization. We also find that the spontaneous spin-polarization survive 
the curvature effect in the BN nanotubes Thus the various nanostructures originating 
from graphitic-BN promise to act as molecular magnets when doped with carbon. Contrast 
to the models proposed by Kusakabe jj| and Choi j^j in which the spin polarization only 
occurs for certain particular structural configurations, the spin polarization in graphitic-BN 
is induced by substitution with arbitrary atomic configuration, indicating easy accessibility 
of experimental synthesis. 

Boron nitride has important advantages over carbon nanostructures. It is far more re- 
sistant to oxidation than carbon and therefore more suitable for high-temperature applica- 
tions in which carbon nanostructures would burn. Unlike carbon nanostructures, BN based 
nanostructures are insulators, with predictable electronic properties. In addition, magnetic 
nanostructures are of scientifically interesting and technologically important, with many 
present and future applications in permanent magnetism, magnetic recording and spintron- 
ics. The carbon substituted BN nanostructure, with conduction electrons that are 100% 
spin polarized due to the gap at the Fermi level in one spin channel, and a finite density of 
states for the other spin channel, can be an ideal half metallic material and can be useful for 
spintronic applications, such as tunneling magnetoresistance and giant magnetoresistance 
elements. 

In conclusion, we have performed first principle projector augmented wave (PAW) po- 
tential calculation with plane wave basis sets to study the effects of carbon substitution on 
graphitic-BN . Our results show that single carbon atom substitution for any atom results in 
a polarized flat band and a sharp peak in the majority density of state below the Fermi level. 
The spontaneous spin-polarization is independent of site of substitution. Compared to other 
metal-free magnets previously proposed, the carbon induced magnetization in graphitic-BN 
is experimentally accessible and the system can be potentially very useful. 
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FIG. 1: The calculated band structures of the graphitic-RN with a boron atom (a) and a nitrogen 
atom (b) substituted by a carbon atom. The solid lines represent the spin-up branches while the 
dot lines represent the spin-down branches. Fermi level is denoted by the light dot line. 
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Fig 2 

FIG. 2: The magnetization density (spin-up minus spin-down) distribution of the graphitic-UN 
with a boron atom (a) and a nitrogen atom (b) substituted by a carbon atom through the sheet 
plane. The unit: spins/A 3 
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FIG. 3: The magnetization dens^tj 
with a boron atom (a) and a nitr<gg§] 
sheet plane. The unit: spins/ A 3 
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FIG. 4: The charge distribution of the graphitic-UN with a boron atom (a), and a nitrogen atom 
(b) substituted by a carbon atom. The unit: electrons/A 3 
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